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Abstract: We extend a previous phenomenological analysis of photon lensing in an ex¬ 
ternal gravitational background to the case of a massless neutrino, and propose a method 
to incorporate radiative effects in the classical lens equations of neutrinos and photons. 
The study is performed for a Schwarzschild metric, generated by a point-like source, and 
expanded in the Newtonian potential at hrst order. We use a semiclassical approach, 
where the perturbative corrections to neutrino scattering, evaluated at one-loop in the 
Standard Model, are compared with the Einstein formula for the deflection using an impact 
parameter formulation. For this purpose, we use the renormalized expression of the gravi¬ 
ton/fermion/fermion vertex presented in previous studies. We show the agreement between 
the classical and the semiclassical formulations, for values of the impact parameter hy of 
the neutrinos of the order of by ~ 20, measured in units of the Schwarzschild radius. The 
analysis is then extended with the inclusion of the post Newtonian corrections in the exter¬ 
nal gravity held, showing that this extension hnds application in the case of the scattering 
of a neutrino/photon off a primordial black hole. The energy dependence of the dehection, 
generated by the quantum corrections, is then combined with the standard formulation of 
the classical lens equations. We illustrate our approach by detailed numerical studies, using 
as a reference both the thin lens and the Virbhadra-Ellis lens. 
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1 Introduction 

According to classical general relativity (GR) massless particles follow null spacetime geodesics 
which bend signihcantly in the presence of very massive sources. The gravitational lensing 
enforced on their spatial trajectories provides important information on the underlying dis¬ 
tributions of matter and, possibly, of dark matter, which act as sources of the gravitational 
held. 

Several newly planned weak lensing experiments such as the Dark Energy Survey (DES) 
[1], the Large Synoptic Survey Telescope (LSST)[2], both ground based, or from space with 
the Wide-Field Infrared Survey Telescope (WFIRST) [3] and Euclid [4], are expected to 
push forward, in the near future, the boundaries of our knowledge in cosmology. 

In the analysis of the dehection by a single compact and spherically symmetric source, one 
signihcant variable, beside the mass of the source, is the impact parameter of the incoming 
particle beam, measured respect to the center of the source, which determines the size of 
the dehection. It is very convenient to measure the impact parameter (6), which is typical 
of a given collision, in units of the Schwarzschild radius = 2GM, denoted as bh = b/vg- 
In the Newtonian approximation for the external background, this allows to scale out the 
entire mass dependence of the lensing event. 

For an impact parameter of the beam of the order of 10® — 10®, the corresponding dehection 
is rather weak, of the order of 1-2 arcseconds, as in the case of a photon skimming the sun. 
Stronger lensing effects are predicted as the particle beam nears a black hole, with dehec- 
tions which may reach 30 arcseconds or more. These are obtained for impact parameters 
bh of the order of 2 x 10“^. Even larger dehections, of 1 to 2 degrees or a signihcant fraction 
of them, are generated in scatterings which proceed closer to the event horizon [5]. In fact, 
as we are going to show, for closer encounters, with the beam located between 20 and 100 
bh, such angular dehections are around 10“^ radians in size, as predicted by classical GR. 

A high energy cosmic ray of 10-100 GeV will then interact with the held of the source by 
exchanging momenta far above the MeV region, and will necessarily be sensitive to radiative 
effects, such as those due to the electroweak corrections. 

Interactions with such momentum exchanges cannot be handled by an effective Newtonian 
potential, as derived, for instance, from the (loop corrected) scattering amplitude. We re¬ 
call that, in general, in the derivation of such a potential, one has to take into account 
only non-analytic terms in the momentum transfer q. These are obtained from a given 
amplitude and/or gravitational form factor of the incoming particle after an expansion at 
small momentum. The analytic terms in the expansion correspond to contact interactions 
which are omitted from the hnal form of the potential, being them proportional to Dirac 
delta functions. 

As one can easily check by a direct analysis, non-analytic contributions originate from 
massless exchanges in the loops, which approximate the full momentum dependence of the 
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radiative corrections only for momentum transfers far below the MeV region. Therefore, 
the validity of the method requires that the typical impact parameter of the beam, for a 
particle with the energy of few GeV’s, be of the order of 10® Schwarzschild radii and not less. 
For such a reason, if we intend to study a lensing event characterized by a close encounter 
between a cosmic ray and a black hole, we need to resort to an alternative approach, which 
does not suffer from these limitations. 

Finally, with the photon sphere located at bh ~ 2.5 for a Schwarzschild metric, one expects 
that very strong deflections are experienced by a beam for scattering events running close 
to such a value of the impact parameter. This is also the radial distance from the black 
hole center at which the scattering angle diverges. A simple expansion of the Einstein 
formula for the deflection shows that this singularity is logarithmic [5]. In such extreme 
cases the beam circulates around the source one or more times before escaping to infinity, 
generating a set of relativistic images [6]. This is also the region where the simple New¬ 
tonian approach, discussed in [5], fails to reproduce the classical GR prediction, as expected. 


1.1 Comparing classical and semiclassical effects 

The analysis of possible extensions of the classical GR prediction for lensing, with the 
inclusion also of quantum effects in the interaction between the particle source and the 
deflector (lens), has not drawn much attention in the past, except for a couple of very original 
proposals [7, 8]. While these effects are expected to be small, even for huge gravitational 
sources such as massive/supermassive black holes, they could provide, in principle, a way 
to test the impact of quantum gravity and of other radiative corrections to the propagation 
of cosmic rays. Glose encounters of a beam with a localized source, which could be a large 
black hole or a neutron star, are expected to be quite common in our universe, although the 
probability of identifying a lensing event characterized by a close alignment between the 
source, the lens and an earth based detector, especially for neutrinos, is exceedingly rare 

[9] . The situation might be more promising for photons in close encounters with primordial 
black holes, revealed by resorting to spaceborne detectors. 

Such is the FERMI satellite [11], with source beams given by Gamma Ray Bursts (GRBs) 

[10] , which could detect fringes between primary and secondary paths of the GRBs on its 
ultra sensitive camera, generated by a gravitational time delay. This approach was termed 
in [10] "femtolensing", due to the size of the Einstein radius characteristic of these events, 
which was estimated to be of the order of a femtoarcsecond. As shown in [10], a classical 
GR analysis based on the thin lens equation can be applied quite straightforwardly also to 
this extreme situation. 

An important point which needs to be addressed, in this case, concerns the quantum features 
of these types of lensing events, since the Schwarzschild radius of a primordial black hole, 
for a gamma ray photon, is comparable to its wavelength. Our analysis draws a path in 
this direction. 

The classical deflections of photons, as pointed out in the past and in a recent work 
[5], can be compared at classical and quantum levels by equating the classical gravitational 
cross section, written in terms of the impact parameter of the incoming photon beam, to 
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the perturbative cross section. The latter is expanded in ordinary perturbation theory with 
the inclusion of the corresponding radiative corrections. The result is a differential equation 
for the impact parameter of the beam, whose solution provides the link between the two 
descriptions. In particular, the energy dependence, naturally present in the cross section 
starting at one-loop order, allows to derive a new formula which relates to the energy E 
of the beam and to the angle of deflection a, bh{E, a). This dependence, which is absent in 
Einstein’s formula, propagates into all the equations for the usual observables of any lensing 
process: magnifications, cosmic shears, the light curve of microlensing events and Shapiro 
time delays. Clearly, such a dependence implies, as noted in [19], that radiative corrections 
induce a violation of the classical equivalence principle in General Relativity. The violation 
of the equivalence principle, viewed from a quantum perspective, is not surprising, since this 
principle is inherently classical and requires the localization of the point particle trajectory 
on a geodesic. It can be summarized in the statement that an experiment will not be able 
to determine the nature of the point particle which is subjected to gravity, except for its 
mass. The notion of a point particle clearly clashes with the quantum description, which 
is, on the other hand, inherently tight to Heisenberg’s indetermination principle. For this 
reason, one expects that the inclusion of radiative corrections will cause a violation of such 
principle. 

Gravity, in this approach, is treated as an external background and the transition ampli¬ 
tude involves on the quantum side, in the photon case, the TVV vertex, where T denotes 
the energy momentum tensor (EMT) of the Standard Model and V the electromagnetic 
current. In the fermion case (f), the corresponding vertex is the Tff, with / denoting 
a neutrino. The comparison between the classical and the semiclassical formula for the 
deflection derived by this method can then be performed at numerical level, as shown in [5] 
for the photons. The energy dependence of the bending angle, for a given impact parameter 
of the photon beam, though small, is found to become more pronounced at higher energies, 
due to the logarithmic growth of the electroweak corrections with the energy. 

The goal of our present work is to propose a procedure which allows to include these 
effects in the ordinary lens equations, illustrating in some detail how this approach can be 
implemented in a complete numerical study. We mention that our semiclassical analysis is 
quite general, and applies both to macroscopic and to microscopic black holes. In the case of 
macroscopic black holes the procedure has to stop at Newtonian level in the external field. 
In fact, post-Newtonian corrections, though calculable, render the perturbative expansion 
in the external (classical) gravitational potential divergent, due to the macroscopic value 
of the Schwarzschild radius. On the other hand, in the case of primordial black holes, the 
very same corrections play a significant role in the deflection of a cosmic ray, and bring to 
a substantial modification of the classical formulas. 

1.2 Organization of this work 

In the first part of our work we will extend a previous analysis of photon lensing [5], 
developed along similar lines, to the neutrino case, presenting a numerical study of the 
complete one-loop corrections derived from the electroweak theory. The formalism uses a 
retarded graviton propagator with the effects of back reaction of the scattered beam on 
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the source not included, as in a typical scattering problem by a static external potential. 
In this case, however, because of the presence of a horizon, we search for a lower bound 
on the size of the impact parameter of the collision where the classical GR prediction and 
the quantum one overlap. Indeed, above the bound the two descriptions are in complete 
agreement. As already mentioned above, both in the fermion as in the photon case [5], this 
bound can be reasonably taken to lay around 20 b^, which is quite close to the horizon of 
the classical source. 

For smaller values of bh (4 < < 20), the two approaches are in disagreement, since the 

logarithmic singularity in the angle of deflection, once the beam gets close to the photon 
sphere, starts playing a significant role. This is expected, given the assumption of weak 
field for the gravitational coupling, which corresponds to the Newtonian approximation in 
the metric. 

The second part of our work deals with the implementation of the semiclassical deflection 
within the formalism of the classical lens equations. We use the energy dependence of 
the angular deflection to derive new lens equations, which are investigated numerically. 
We quantify the impact of these effects both in the thin lens approximation, where the 
trigonometric relations in the lens geometry are expanded to first order, and for a lens with 
deflection terms of higher order included. As an example, in this second case, we have 
chosen the Virbhadra-Ellis [12] lens equation. The observables that we discuss are limited 
to solutions of these equations and to their magnifications, although time delays, shears and 
the light curves of a typical microlensing event can be easily included in this framework. We 
anticipate that the effects that we quantify are small and cover the milliarcsecond region, 
remaining quite challenging to detect at experimental level. We hope though, that the 
framework that we propose can draw further interest on this topic in the future, both at 
theoretical and at phenomenological level. 

In the third part of our study we discuss the post Newtonian formulation of the impact 
parameter formalism, and apply it to the case of a compact source with a microscopic 
Schwarzschild radius. This is the only case in which the gravitational corrections to the 
Newtonian cross section can be consistently included in our approach in a meaningful 
way. We then summarize our analysis and discuss in the conclusions some possible future 
directions of possible extensions of our work. 

2 Gravitational interaction of neutrinos 

We start our analysis with a brief discussion of the structure of the gravitational interaction 
of neutrinos, building on the results of [13, 14], to which we refer for additional details, and 
that we are going to specialize to the case of a massless neutrino. An analysis of gravity 
with the fermion sector is contained in [15]. We simply recall that the dynamics of the 
Standard Model in external gravity is described by the Lagrangian 

d^Xy/^CsM + X J RH"'!!. (2. 

This includes the Einstein term Sq, the Ssm action and a term Sj involving the Higgs 
doublet H [16], called the term of improvement. Ssm, instead, is obtained by extending 


S = Sg + Ssm + Si = -^ d^x^R + 
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the ordinary Lagrangian of the Standard Model to a curved metric background. The term 
X is a parameter which, at this stage, is arbitrary and that at a special value (x = Xc = 1 /6) 
guarantees the renormalizability of the model at leading order in the expansion in k. 
Deviations from the flat metric will be parametrized in terms of the 

gravitational coupling k, with k? = IGttG and with G being the gravitational Newton’s 
constant. At this order the metric is given as + nh^iy, with describing its 

fluctuations. We will consider two spherically symmetric and static cases, corresponding to 
the Schwarzschild and Reissner-Nordstrom metrics. The first, in the weak field limit and 
in the isotropic form is given by 


ds^ 


1 - 


2GM\ 


dr - [1 + 


2GM\ 


In this case the fluctuation tensor takes the form 




2GM - 

I ^1 


dx ■ dx. 


S/iu — 'Hfiv 


K\X\ 


flUl 


( 2 . 2 ) 


(2.3) 


The inclusion of higher order terms in the weak field expansion will be discussed in the 
following sections. 

The coupling of the gravitational fluctuations to the fields of the Standard Model involves 
the EMT, which is defined as 


^ _ 2 d{SsM + Si) 

with a tree-level coupling summarized by the action 

Sint = -^ j d^xTg.hf^^', 


(2.4) 


(2.5) 


where is symmetric and covariantly conserved. The complete expression of the EMT of 
the Standard Model, including ghost and gauge-fixing contributions can be found in [17]. 
The Higgs field is parameterized in the form 


H = ( I (2.6) 

in terms of h, cj) and 0^, which denote the physical Higgs and the Goldstone bosons of 
the Z and W 's respectively, v is the Higgs vacuum expectation value. The terms of the 
Lagrangian 5/, generate an extra contribution to the EMT which is given by 


= -2x(9'^9^ - = -2x{d^d'^ - r/^^D) (y + y + /i) , (2.7) 


the term of improvement, which can be multiplied by an arbitrary constant (x). As men¬ 
tioned above, it is mandatory to choose the value x = 1/6 for any insertion of the EMT on 
the correlators of the Standard Model. These are found to be ultraviolet finite only if 
is included [16-18]. 


- 6 - 







We will be dealing with the Tff vertex, where T denotes the EMT and f = Uf a neutrino 
of flavour /, and work in the limit of zero mass of the neutrinos. The vertex, to lowest 
order, is obtained from the EMT of the neutrino. For instance, the explicit expression of 
the EMT for the (left-handed, v = ul) electron neutrino is given by 


with 


- 2i 


2 e 


1 - 7 ' 




\/2 sin Ow 

I/) I - , 


sin26lw^^^ 2 

W+ + W-^ 


1-Y 




+ 


sin 2'(9w 

^ eW^ + 67 ^ ^ ^ W~ ). 


^/2 sin y ' 2 ^ '2 

In momentum space, in the case of a massless fermion, the vertex takes the form 


( 2 . 8 ) 


(2.9) 


^ +p^y + ^-(p^ +p2)m _ +^^)). (2.10) 

while in the case of neutrinos we have 


YiO)llU ^yiO)l,U ( 2 . 11 ) 

with Pl = (1 — 75)72 being the chiral projector. We refer to appendix G for a list of the 
relevant Feynman rules necessary for the computation. 

We will denote with 

fm^lu ^ u{p^)V^^)^^''u{pl), ( 2 . 12 ) 

the corresponding invariant amplitude, a notation that we will use also at one-loop level 
in the electroweak expansion. We introduce the two linear combinations of momenta p = 
Pl + P 2 and q = Pl — P 2 express our results. It has been shown that the general Tff 
vertex, for any fermion / of the Standard Model, decomposes into six different contributions 
[13], but in the case of a massless neutrino only three amplitudes at one-loop level are left, 
denoted as 


^ rjnflh' ^ rj^fJ,U 


iCT- 


(2.13) 


In the expression above, the subscripts indicate the contributions mediated by virtual Z 
and W gauge bosons, while CT indicates the contribution from the counterterm. 

We show in Fig. 1 some of the typical topologies appearing in their perturbative ex¬ 
pansion. 

Two of them are characterized by a typical triangle topology, while the others denote terms 
where the insertion of the EMT and of the fermion field occur on the same point. The 
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Figure 1. The one-loop Feynman diagrams of the neutrino vertex in a gravitational background. 
The dashed lines can be Z and W. 


computation of these diagrams is rather involved and has been performed in dimensional 
regularization using the on-shell renormalization scheme. Neutrinos interactions, in the 
limit of massless neutrinos, involve only few of the structures of the Tff tensor decompo¬ 
sition presented in [13]. In this case we are left with only one tensor structure and hence 
only one form factor for each sector 

Tw = i {Q^,mf,mw)HP2)0^\u{pi), (2.14) 

where we have defined the vertex 


0^1 = 


(2.15) 


The counterterms needed for the renormalization of the vertex can be obtained by promoting 
the counterterm Lagrangian of the Standard Model from a flat spacetime to the curved 
background, and then extracting the corresponding Feynman rules, as for the bare one. We 
obtain 




where we have denoted with the neutrino self-energy 


S^(p^) = 


Gf 


P'zip^) + 'P‘wip 


L 


167r2y/2 

which is a combination of the self-energy contributions 


(2.16) 


(2.17) 


S^(p 2 ) = -4 

'P‘W) = - 2 m 


(m2 -p 2m%/) Bi {p^,m 


2 

z 


2 Bi [p^, 0, m^) -|- 1 


2 

/’ 


“I" 


(2.18) 

(2.19) 


with 


Pi {p‘^,ml,ml) = 


mt — m, 


2 r 


2 p 2 


Boip^, ml, mj) - ^ 0 ( 0 , ml, mj) 


- 2^o(/,mg,m?),(2.20) 
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expressed in terms of the scalar form factor Bq, given in appendix E together with all the 
other relevant scalar integrals. We have denoted with mz and mw the masses of the Z and 
W gauge bosons; with the virtuality of the incoming momentum of the EMT and is 
the mass of the fermion of flavor / running in the loops. 

The explicit expressions of the form factors appearing in (2.14) is given by 


7 9 4 Tfi ry 

h = -2mz-^ + 


2 + 


7 m\ 
2>q^ 


Ao{m%) - 


17 7 m\ 4m| 

6 


i3o((?^0,0) 


+ i "^1(2 + q^) (3m| + 2q^) ^o(g^ m|, m|) 

6 q 

+ 9^) ^o(0, m|) — rn^z + 9^)0,0(?.21) 


with Co denoting the scalar 3-point function, and with the form factor related to the 
exchange of the W s given by 


/r = -4- - (m) + -2m^) - ^{m) + 2rn^) {Ao{mj) - 

- {m^r + q'^4 BoiO,mj,m‘^) + “ 24m® - lOm^g^ 

-|- mj (72 m^ -|- 46 m^ q'^ + q^) — 2 rn^ (24 m^ -|- 42 m^ (7^ -|- 17 5^)^ BQ{q'^, ^/) 

+ 2 ^4 ^12 m® -|- 12 mj- -|- 4 m^ (2 m^ -|- 9^)(3 m^ -|- 2 <7^) 

+ (—36 m^ — 16m^ ^0(9^) '^w) 2^ (^/ “ (^/ 2 

^ (3m^ -4mjm^-F5mj m^-h 4 m^) ^ Co(m|, m^, m^) 

+ ^ ^4 m® -|- m® (g^ — 4 rn^) — 2 m^ “I" 9^)^ (4 + <7^) 

- m) (2 m^ (6 m^ -h Co(m^, rrij, mj) 

9 

TTh^ 

+ (20m^ -|- 25 m|y (7^ -|- 6m^ (7^) Co(m^, m^, (2.22) 


Being the computations rather involved, the correctness of the results above has been 
secured by appropriate Ward identities, whose general structure has been discussed in [17]. 
As an example, by requiring the invariance of the generating functional of the theory under 
a diffeomorphic change of the spacetime metric, one derives the following Ward identity 






where is the fermion two-point function, diagonal in flavor space [13]. From this 

equation one obtains 


0 = /f - ^S|( 0 ) 
o = /r-^s^(o). 
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(2.24) 







Soul'ce 


Observer 



Figure 2. The deflection of the trajectory of a massless particle P approaching a black hole. 


which, as one can check, are identically satished by the explicit expressions of fz and fw 
given above. 

In the case of MeV neutrinos, the expressions of the two form factors simplify considerably, 
since the typical momentum transfer slv?{6/2) may be small. These expansions, 

in fact, are useful in the case of scattering and lensing of neutrinos far from the region of the 
event horizon, of the order of 10^ — 10® horizon units. As we are going to see, an expansion in 
provides approximate analytical expressions of the hh{a) relation, connecting the impact 
parameter to the angle of deflection a, valid at momentum transfers which are smaller 
compared to the electroweak scale, i.e. q^ <C 1. We will come back to illustrate this 
point more closely in the following sections. 

In these cases the expression of the renormalized fz form factor takes the form 


pZ (ren) _ 

low q 

while the W form factor is slightly lengthier 


11 

18 


q , 


(2.25) 


nW (ren) _ 

low q 


5 — 98 + 243 — 194 mjm^ + 44 


4.^4 




36 {m'j — rri^Y 

+6 ^lOm^m^ — In 


(2.26) 


3 Cross Sections for photons, massive fermions and scalars 

Before coming to a discussion of the 1-loop effects in the scattering of neutrinos, we briefly 
summarize the result for the leading order cross sections for fermions, photons and scalars 
using in an external static background [19] [13, 14]. We just recall that the scattering matrix 
element is written as 

iSif = -^ y d^x{p 2 \hf,u{x)T>"''{x)\pi), (3.1) 

where V is the integration volume where the scattering occurs, which gives 

{p2\h^.{x)T>^''{x)\p,) = Ku{x)i^{P2)V^''i’{pi)e^^-^. (3.2) 

Denoting with i and / the initial and hnal neutrino, we have introduced plane waves 
normalized as 

'Ipi{pi) = Niu{pi), = u{pi)u{pi) = l, (3.3) 
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and similarly for 'ipf, while V denotes a finite volume. The Ei {E 2 ) are the energy of the 
incoming (outgoing) particle respectively. 

In momentum space the matrix element is given by 

tSf, = (3.4) 

in terms of the gravitational fluctuations in momentum space h^i,{q). For a static external 
held the energies of the incoming/outgoing fermions are conserved {Ei = E 2 ^ E). 

The Fourier transform of in momentum space is given by 


hf,u{qo,^ = j 

which for a static held can be expressed as 


(3.5) 


{qo,q) = 27r(5(go)Vi.(€), 

in terms of a single form factor hQ{q) 


hiMui^ = ho{q)S, 


fiiy 


with ho{q) = 


kM 


2q 




The squared matrix element in each case takes the general form 

= i6v^^27r<5(go) r 2 :^^"""(P 1 ,P 2 ) V(9) V(?) 


(3.6) 


(3.7) 


(3.8) 


where T is the transition time. Specihcally, in the case of a massive (Dirac) fermion one 
obtains 

= tr [( 1^2 + {pi,p2){i>i + m)V^{pi,p2)] , (3.9) 

where the vertex is in this case given by 


VmiPi,P 2 ) = ^(j^iPi +P 2 y + l’'{pi +P 2 T - - 2m)) 

which gives a cross section 

2 


da 


( 0 ) 

/ 


sm 


GM Y 


9 „ , 1 m? 1 3 9 „ , 

cos^ i?/2+ vTvrry + 77^771 + 7 77777 cos^ 19/2 

4|pi|2 4 |pi|4 4|pi|2 


(3.10) 


(3.11) 


In the case of a neutrino, the corresponding cross section is obtained by sending the fermion 
mass m of the related Dirac cross section to zero, giving 


da 

dQ 


( 0 ) 


GM 


2 ^ 
cos 2 , 


(3.12) 


which is energy independent. Notice that the inclusion of the chiral projector Pl in the 
expression of the neutrino amplitude, which carries a factor 1 / 2 , makes the neutrino and 
Dirac cross sections coincide. The same 1/2 factor, in the Dirac case, appears in the average 
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over the two states of helicity, while the axial-vector terms induced by Pl are trivially zero 
(see [20] for typical studies of polarized processes). 

Ill the photon case one obtains 

Ai,A 2 

(3.13) 

where denotes the polarization vector of the photon, with an interaction vertex which is 
given by 

k2) = i| (All • k2) fc2)| , (3.14) 

where 


C^i/pa — ^p<^ ? 


Dpi/pa {h\ , A^2) — Vpu ^1 (T ^2 p 


+ r}^pkiak2u - Vpa kipk2u + (fJ-^ i^) 


The cross section for a photon is then given by 


da 

dO, 


( 0 ) 


= {GMycoP{e/2). 


(3.15) 


Finally, in the case of a scalar the relative expression is given by 

J’"^^'"(pi,P2) = VY{pi,p2)VP^{pi,p2) , 


(3.16) 


with 


vr = -i {pipP2.C^^^‘^ - 2x [(pi +P2)^ (Pi +P2)" - p'^"(pi +P2)'] } , (3.17) 


where we have included the minimal and the term of improvement [17]. For a conformally 
coupled scalar y = 1/6. The cross sections, in this case, are given by 


da _ f (GM)2 csc^( 0 / 2 ) x = 0 

~ i (T)'cot4(0/2) x = l/6 


We show in Fig. 3 the expressions of these three cross sections at different energies, normal¬ 
ized by l/{2GMy and denoted as a. In panel (a) we consider the scattering of a massive 
fermion, together with the massless limit, which applies in the neutrino case. We have in¬ 
cluded in (b) and (c) two enlargements of (a) which show how the massive and the massless 
cross sections tend to overlap for energies of the order of 1 GeV. In panel (d) we show the 
cross sections for the photon (s = 1), for the neutrino (s = 1/2) and for the conformally 
coupled scalar (s = 0). 
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dtr/dfi 
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Figure 3. Normalized (ct = cr/(2GM)^) cross sections for massive and massless fermions. In the 
massive case m is the electron mass (a). Two enlargements of (a) are in (b) and (c). Panel (d) 
shows the cross sections for photons (s = 1), massless neutrinos (s = 1/2) and conformally coupled 
scalars (s = 0). 


3.1 The neutrino cross section at 1-loop 


In the neutrino case, at 1-loop level, Eq. (3.12) is modified in the form 

4Gf 


— = c23^2 Cos^^/2 
dO. sin^ 0/2 


1 + 


16 7r2y/2 . 




(3.19) 


whose explicit expression has been given in appendix C. In the massless approximation for 
the neutrino masses, loop corrections do not induce flavor transition vertices, such as those 
computed in [14]. 

In the case of neutrinos of an energy E in the MeV range, the expression above simplifies 
considerably and takes the form 


da 

dQ 


= G^ 


^M" C°s)»/2 ji + _gF 


sin"^ 6/2 


vr" 


V2 


11 


+ 


5 — 98 + 243 


4^4 


18 ' 36 {m'j — m^)^ 


— 194 mjm^ + 44 + 6 ^10 — 15 rrij-m'^ + 2 In 


m. 


m 


w 


E'^ sin^ 


(3.20) 


We show in Fig. 4 three plots of the tree level and one-loop cross sections for an energy of 
the incoming neutrino beam of 1 MeV, for 2 different angular regions (plots (a) and ( 6 )), 
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Figure 4. Differential cross section for MeV neutrinos in units of with Ts the Schwarzschild 
radius. 


together with a global plot of the entire cross section (plot (c)) for the rescaled differential 
cross section da/dQ = 1/r^ da/dQ. Notice that the tree-level and one-loop results are 
superimposed. We can resolve the differences between the two by zooming-in in some 
specific angular regions of the two results, varying the energy of the incoming beam. The 
result of this analysis is shown in Fig. 5, where in plots (a) and (6) we show the rescaled 
cross section da/dO. as a function of the scattering angle 0, for three values of the incoming 
neutrino beam equal to 1 GeV, 1 TeV and 1 PeV. PeV neutrinos events are rare, due to 
the almost structureless cosmic ray spectrum, which falls dramatically with energy. They 
could be produced, though, as secondaries from the decays of primary protons of energy 
around the GZK [21, 22] cutoff, and as such they are part of our analysis, which we try to 
keep as general as possible. 

It is clear from these two plots that the tree-level and the one-loop result are superimposed 
at low energies, with a difference which becomes slightly more remarked at higher energies. 
A similar behaviour is noticed in the cross section for scatterings at larger angles. Also 
in this case the radiative corrections tend to raise as the energy of the incoming beam 
increases. This behaviour is expected to affect the size of the angle of deflection a as we 
approach the singular region of a black hole. In fact, a is obtained by integrating the 
semiclassical equation (4.1), introduced below, and large deviations are expected as the 
impact parameter reaches the photon sphere. As we are going to illustrate in the next 
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Figure 5. Differential cross section: tree level and one-loop contribution for a wide range of 
energies. 


sections, the bh{E, a) relation is significantly affected by the behaviour of the cross section 
at large 9 as bh ^ 3/2rs. This is the closest radial distance allowed to a particle approaching 
the black hole from infinite distance without being trapped. Therefore, these differences in 
a for large 9 are going to render bh sensitive on the changes in energy of the neutrino beam 
for such close encounters of the neutrinos with a black hole. 


4 Impact parameter formulation of the semiclassical scattering 

As pointed out in previous studies [5, 7, 14, 23], the computation of the angle of deflection 
for a fermion or a photon involves a simple semiclassical analysis, in which one introduces 
the impact parameter representation of the specific classical cross section and equates it to 
the quantum one. The classical/semiclassical scattering process is illustrated in Fig. 2, with 
a denoting the angle of deflection. By assuming that the incoming particle is moving along 
the z direction, with the source localized at the origin, and denoting with 9 the azimuthal 
scattering angle present in the quantum cross section, we have the relation 


b 

db 

da 

sin0 

d9 

dQ 


between the impact parameter b and 9, as measured from the z-direction. This semiclassical 
equation [7, 23] allows to relate the quantum and the classical features of the interaction 
between the particle beam and the gravitational source. The explicit expression of b{a), at 
least for small deflection angles, which correspond to large values of the impact parameter, 
can be found either analytically, such as at Born level and, for small momentum transfers 
also at one-loop, but it has to be obtained numerically otherwise. The solution of (4.1) 
takes the general form 

bl{a) = bl{9) + 2 J d9'sm9'-^, (4.2) 

with bf^{9) denoting the constant of integration. The semiclassical scattering angle a is 
obtained from (4.2) as a boundary value of the integral in 9 of the quantum cross section. 
As discussed in [5], the integration constant derived from (4.2) has to be set to zero (for 
9 = it) in order for the solution of (4.1) to match the classical GR result for a very large b^- 
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Figure 6. Plot of a, the angle of deflection, versus the impact parameter bh for 20 < bh < 100, 
having inverted the b{a) solution of (4.1) at various 1/6" orders. 


In the case of a point-like gravitational source and of neutrino deflection, one obtains 
from (4.2) the differential equation 

2 


le 


= -2 


GM 
,2 e 


cos^ - sin 0 . 


(4.3) 


Notice that the variation of b with the scattering angle 9 is negative, since the impact 
parameter decreases as 9 grows, as we approach the center of the massive source. A com¬ 
parison of this expression with the analogous relation in the photon case ( 7 ) shows that the 
two equations differ by a simple prefactor 


~d9 


1 db"^ 


cos^ 2 d9 


db'^ 

-M 


Q 

= — 2 G^M^ cot'^ - sin 0 . 


The solution of (4.3) takes the form 

b^{a) = 4G^M^ ^ ^ ^)) ’ 

and in the small a. (i.e. large b) limit takes the asymptotic form 


b = GM ( - + -(l + ln8-31na) ) +0(a^) 
a 3 


(4.4) 


(4.5) 


(4.6) 


which allows us to identify the deflection angle as 

GM 
a ~ 4—;— 


(4.7) 

in agreement with Einstein’s prediction for the angular deflection. This is the result ex¬ 
pected from the classical (GR) analysis. The inversion of the asymptotic expansion (4.6) 
generates the asymptotic behaviour 


« = ^ + h) + 4(1^' + Oil/bl) 


bh bl 


bl 


(4.8) 
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Massless Particles 




(b) 

Figure 7. Angle of deflection as a function of the impact parameter: (a) spin dependence; (b) 
mass dependence. 


which corresponds to the general functional form 

« = 7—b 2fc+i (®2fc+i + di In bh + d2 In^ bh + ■ ■ ■ + dk In^ bh) ■ (4-9) 

k>i k>l 

As shown in Fig. 6, the analytic inversion of (4.6), given by (4.8), is very stable under an 
increase of the order of the asymptotic expansion over a pretty large interval of b^, from 
low to very high values. Solutions (4.8) and (4.9) can be obtained by an iterative (hxed 
point) procedure, which generates a sequence of approximations ao —t ai —t ...—)• to 
a{bh) implemented after a Laurent expansion of (4.6) and the use of the initial condition 
ao = 2/6/i. The approach can be implemented also at one-loop and with the inclusion of 
the post-Newtonian corrections, if necessary. 

The logarithmic corrections present in (4.9) are a genuine result of the quantum ap¬ 
proach and, as we are going to discuss below, are not present in the classical formula for the 
deflection. Radiative and post-Newtonian effects, not included in (4.8), give an expression 
for a{bh) which coincides with the form (4.9), with specihc coefficients {an,dn) which are 
energy dependent. This is at the origin of the phenomenon of light dispersion (gravitational 
rainbow) induced by the quantum corrections, which is absent at classical level [19]. 
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Eq. (4.9) will play a key role in our proposal for the inclusion of the radiative corrections 
in the classical lens equation. Such equation will relate the angular position of the source 
in the absence of lensing, /3, to a(b). 

We give, for completeness, the analogous expressions in the case of the scalar and for a 
massive fermion. For a massless scalar we have the relation 


a = 


3 6 /^ 


bl 


12 ln3 — 1 4 1 5 

^ ) + 0 ( 1 / 6 J, 


(4.10) 


while for a massive fermion the corresponding expression becomes more involved and takes 
the form 




1 1 


a = 




-m)){AE^ -2 E^mj + mj^ 
2 


4E^-2 E^mj + mj bh hi 
X (^)+^E^{l+\n^) + E^m)\nU-QE^{AE‘^ +rr?)\n 
AE^iAE^ + mj) 


m% 

1_4 


+ 


8E^ — 8 E^m?j + 4 E'^mjr — 


In bh 


+ o{i/bl), 


(4.11) 


where E and ruf are the energy and the mass of the fermion respectively. One can easily 
check that in the limit E ^ nif Eq. (4.11) reproduce the formula for the massless fermion 
(neutrino). We have plotted the behaviour of the formulas for the deflection in Fig. 7. 
As one can immediately notice from the plots presented, the angular deflection is much 
less enhanced in the scalar case compared to the remaining cases, showing a systematic 
difference respect to the classical prediction form Einstein’s deflection integral (F.l). The 
angular deflection in the scalar case is signihcantly affected by the choice of x the free 
coupling factor of a scalar held to the external curvature R. We have chosen in panel (a) 
the two cases of x = 0 (minimal coupling) and of conformal coupling (x = 1/6) as typical 
examples. 


4.1 Bending at 1-loop 

Moving to the one-loop expression given in (3.19), we can derive an analytic solution of the 
corresponding semiclassical equation (4.1) for b = b{E,a), in the limit of small momentum 
transfers. For this reason we perform an expansion of (3.19) in up to 0((g^/m^)^) 

and solve (4.1) in this approximation for bl{E,a), obtaining 


bl{E,a) 


— 1 -|- CSC — 2 In l^sin 

-F20 D 2 {E) In (^sin - 
—2 D 2 {E) cos 2a In ^sin^ 



+ Ci{E) 


4^2 (if) cos a 

^)-E2iE), 


1 -|- cos a -|- 4 In |^sin — j 
— 8 D 2 {E) cos a In ^sin^ — 


-7 C2{E) cos ^ ^ 

^ — G 2 { E ) cos 2q: 


(4.12) 


with the coefficients C,D, E and G are functions of the energy and of the masses of the 
weak gauge bosons. Their explicit expressions can be found in appendix D. The impact 
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Figure 8. Plots of the impact parameter bh versus a, the angle of deflection, for 20 < bh < 100 for 
the classical and quantum solution. 

parameter bh{a), as shown in the same appendix, has a dependence on the angular deflection 
a which can be summarized by an expression of the form 

2 

a) = —h c{E) a + d{E) a ln(Q:) + f{E) + g{E) In a + h{E) In^ a + O(a^) 
a 

(4.13) 


that we can invert in order to get a(E, bh)- This is given by 


a{E,bh) 

A{E) 


2 1 

bh bl 


{2 + 4Ci{E))\ogbh + A{E)] +0{l/bl) 


-2 Ci{E) - C2{E) + E2{E) + 4E2iE) + G2{E) + -. 

O 


(4.14) 


We show in Fig. 8 some plots of the impact parameter bh as a function of the deflection 
angle in a range closer to the horizon of a black hole, computed using the Newtonian ap¬ 
proximation derived from the metric (2.2). The region involved covers the interval between 
20 and 100 horizons. The numerical results refer to the GR solution and to the full one- 
loop prediction respectively. The classical expression and the quantum one start differing 
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as we approach the value of 6 /^ ~ 20 , and are characterized by a certain dependence on the 
energy of the incoming beam. Shown are the plots corresponding to neutrinos of energies 
in the TeV and the PeV range respectively. In these regions the lensing is very strong, cor¬ 
responding to 10^ arcseconds and larger. As the neutrino (or the photon) beam gets closer 
to the photon sphere (xq = 3/2rs), which is the point of maximum approach, the angular 
deflection diverges. This is the impact parameter region where one expects the formation of 
relativistic images. The divergence can be parameterized by an integer n, with = Inn, 
and n tending to infinity. The integer is the winding number of the beam path around the 
photon sphere. In the external neighborhood of the point of closest approach the beam still 
escapes to infinity, forming an infinite set of images which are parameterized by the same 
integer n [ 6 ]. 


5 1/b'^ contributions to the deflection 


It is interesting to compare the classical GR prediction for the deflection with the result of 
(4.9), by resorting to a similar expansion for the deflection integral. This has been studied 
quite carefully in the literature, especially in the limit of strong lensing [24, 25]. The 
1 /b'^ expansion has been shown to appear quite naturally in the post-Newtonian approach 
applied to the Einstein integral for light deflection. 

We recall that Einstein’s expression in GR is given by the integral 


a(ro) 




and can be re-expressed in the form 


(5.1) 


2 f \ . <^1 . 

Jo y/l- 2s - y‘^ + 2sy^ 


- TT, 


(5.2) 


with the variable s = rs/(2ro) being related to the ratio between the Schwarzschild radius 
and the distance of closest approach between the particle and the source, tq. The exact 
computation of this integral is discussed in appendix F, and involves elliptic functions. 
Additional information on a(ro) is obtained via an expansion of the integrand in powers of 
s and a subsequent integration. This method shows that the result can be cast in the form 


a{hh) 


^ ^ ^ ^ ^ 
^ ^ ^ ^ 


(5.3) 


with 


ai = 2, 


15 

ao = -TT, 

16 ’ 



04 


3465 

-TT, 

1024 



(5.4) 


The coefficients ai differ from those given in [25] (up to 07 ) just by a normalization. They 
are obtained by re-expressing s = s(ro) in terms of the impact parameter bh using the 
relation 


= x„ (1 - - j (5.5) 
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between the impact parameter and the radial distance of closest approach, having redefined 
xq = vq/{2GM). This can also be brought into the form 


2hh 1 

Xo = cos - cos 
v 3 3 


33/2 \ 

J 


(5.6) 


An expression equivalent to (5.6) can be found in [5]. Eq. (5.6) can be given in a 1/bh 
expansion 


xo = bh 


3 

8bh 


^bl 


105 




(5.7) 


which will turn useful below. 

We can invert (5.3) obtaining the relation 


bh {( y ) =—+ ^2) + TX (®2 ~ 302 03 + 2 04) 

a 2 o '16 

q,3 

H- (8 05 — 16 02 04 — 8 a| + 20 0^03 — 50^) + C>(a^), (5-8) 

128 ^ 

which differs from (4.9) by the absence of logarithmic terms in the impact parameter bh 
and by the energy independence of the coefficients. The inclusion of the extra contributions 
mentioned above, in the classical GR expression, becomes relevant in the case of strong 
lensing. We show in appendix B how the inclusion of the additional 1/6)) terms in the 
expansion of the angular deflection can be extended to the case of a continuous distribution 
of sources/deflectors. This provides a simple generalization of the standard approach to 
classical lensing for such distributions. 


6 Lens equations and l/b^ corrections 

The standard approach to gravitational lensing in GR is based on an equation, derived 
from a geometrical construction, which relates the angular position of the image ( 0 /) to 
that of the source (/3), with an intermediate angular deflection (a) generated on the lens 
plane. In this section we are going to briefly review this construction, which is based on the 
asymptotic expression for the angular deflection (a ~ 2,jbh), and discuss its extension when 
one takes into account more general expansions of a{bh) of the form given by Eq. (5.3). 
The extension that we consider covers the case of a thin lens and concerns only the extra 
l/ 6 )( terms derived from classical GR. The discussion is preliminary to the analysis of the 
next section, where we will consider the inclusion of the radiative effects, parameterized by 
(4.9), into the classical lens equation. 

6.1 The lens geometry 

We show in Fig. 9 the lens geometry in the case of a continuous distributions of sources 
and deflectors. A simplified picture of the geometry, with pointlike source and deflector is 
shown in Fig. 10. We indicate with j3 the oriented angle between the optical axis (OP) 
(taken as the 2 ; axis) and the unlensed direction of the source (OS). Oj denotes the angle 
formed by the visual line of the image (01) with the optical axis. We also denote with 
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Figure 9. Geometric construction of the lens for a continuous distribution of sources. Shown 
are the plane S of the source distribution and the plane of the lens L. The line OI identifies the 
direction at which the observer sees the image after the angular deflection a.. 


Optical Axis 



Source Plane 


Figure 10. The thin lens geometric construction where the source S, the lens V and the observer 
O lie on the same plane. Notice that figure is not to scale, since Dql and Uls are far larger than 
the length oiVR. 

Dql the distance between the observer and the lens plane; with D^s the distance between 
the lens plane and the source plane and with Dqs the distance of the source plane from 
the observer, a is the (oriented) angle of deflection, measured clockwise as all the other 
angles appearing in the geometrical construction. We also introduce the relations, valid for 
Dls^Dql much larger than the size of the lens, typical of a linear lens, 

fJ=P’S = pDos SI = &Dls PI = 9iDos- (6.1) 

The thin lens equation follows from the approximate geometrical relation 

PI = PS + SI i.e. 0 =01-11^. (6.2) 

Pos 
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Figure 11. Geometrical construction for the primary Ip and secondary 1^ images generated by 
the two geodesics of the isotropic emission. Shown are the source S, the lens, represented by the 
dotted circle, the observer O and the primary Ip and secondary 1^ angular positions involved in the 
discussion. 


Denoting with ^ a 2-D vector in the lens plane, it is convenient to introduce two scales r]Q 
and ^0 defined as 


V = Vo y 


i = VR = Cox 


Using the lens equation in the geometric relation 


\PI\ _ Dos 
|V7?| Dol 


we find the relation 


^DlsDol 

y = X — Oi-- --— = X — a 

P>os^o 


with 


% _ Dqs 

Co Dls 


(6.3) 


(6.4) 


uDls Dol 


a = a 


Dos Co 


(6.5) 


which defines the thin lens equation. It is possible to give a simpler expression to the 
equation above if we go back to (6.2) and perform simple manipulations on the angular 
dependence. On the lens plane (Fig. 10) the equation takes the scalar form 


(3 = 9j-a ——, ( 6 . 6 ) 

Dos 

which can be extended to the case of stronger lensing by the inclusion of the contributions 
of the 1/6” corrections in a{b). Use of the Einstein relation a = AGM/b and of the relation 
6 ~ OiDol brings (6.6) into the typical form 


^ 2 _ Dls ^GM 
9i ^ Dos Dol ’ 


/3 = 0 / 


(6.7) 










which dehnes the thin lens approximation, with 9e being the Einstein radius. For a source 
S aligned on the optical axis together with the deflector and the observer O (see Fig. 11) 
- which is dehned by the segment connecting the observer, the lens and the plane of the 
source (with /3 = 0) - the images will form radially at an opening Oj = 9 e and appear as a 
circle perpendicular to the lens plane. For a generic /3, instead, the primary and secondary 
image solutions are given by the well-known expressions 

9i± = ^±\{^^ + A9lY'\ ( 6 . 8 ) 

It is quite straightforward to extend this derivation with the inclusion of the 1/6"" correc¬ 
tions in the a{h) relation and test their effect numerically [25]. This is part of a possible 
improvement of the ordinary (quadratic) thin lens equation which can be investigated more 
generally in conditions of strong lensing. In that case one can also adopt an equation which 
includes deflections of higher orders, as we will discuss in the following sections. For the 
moment we just mention that the inclusion of the higher order 1/6” contributions given by 
(5.3) modihes (6.6) into the form 


(n) 


n n 0% ^9 j, 

■' n>2 J 


with 


/ 9 (”) = r'^n 


Dls 


(6.9) 


( 6 . 10 ) 


DosDl^ 

Another observable that we will investigate numerically is going to be the lens magnihcation. 
For this purpose we recall that light beams are subject to deflections both as a whole but 
also locally, due to their bundle structure. Rays which travel closer to the deflector are 
subject to a stronger deflection compared to those that travel further away. This generates 
a difference in the solid angles under which the source is viewed by the observer in the 
unlensed and in the lensed cases. In the simple case of an axi-symmetric lens the ratio 
between the two solid angles can be dehned in the scalar form 

-1 


= 


9/3 sin/3 
d9i sin 9 1 


( 6 . 11 ) 


In the case of a thin lens (6.5), the analogous expression is given by 

For this lens the analysis simplihes quite drastically. Using the expression of the two images 
9i± given in (6.8) one obtains the simple expression for the primary and secondary images 




-1 


where the Einstein angle is dehned as usual 


9e = JAGM 


Dol 


with 


X = 


Dls 

Dos 


(6.13) 


(6.14) 


24 - 









It is convenient to measure the angular variables in terms of the Einstein angle as 
^ = ^/6e, 0 = 9i/6e, with 

= f ±^1 + ^, (6.15) 

then the total magnification takes a rather simple form 




2 + p^ 


(6.16) 


This equation is commonly used to calculate the light curve in the microlensing case. We 
refer to [26] for a short review on this point. 


6.2 Nonlinear effects in strong deflections 

In conditions of strong lensing, the linear approximations in the trigonometric expressions 
are not accurate enough and one has to turn to a fully nonlinear description of the geometry, 
expressed in terms of the angular variables which are involved. We illustrate this point by 
taking as an example a typical lens equation, which in our case is given by the Virbhadra- 
Ellis construction (VE) [27]. 

Following Fig. 10, we recall that the VE lens equation is based on the geometrical relation 
[27, 28] 


PS = PI- SI, 


(6.17) 


which gives 


tan/3 = Dos^Si,n9i — Zli5(tan0/ + tan(a — 9i)), (6.18) 

under the assumption that the point R in Fig. 10 lies on the vertical plane of the lens. 0/ 
is the angle at which the image is viewed by the observer and /3 is the unlensed angular 
position of the source. Within this approximation we can use the geometric relation 

b = DoLsin9j, (6.19) 

which allows to relate the image position 0/ to the angular deflection of the beam a. Notice 
that this approximate relation is justified by the fact that the distances Pql and ^OS are 
very large compared to the radius of closest approach tq. In this limit the two segments 
VH and VT are treated as equal. 

We remind that (6.18) is not the unique lens equation that one can write down, but, 
differently from Eq. (6.5), it can be used in the case of strong lensing. It takes into account 
the nonlinear contributions to the angular deflection by the introduction of the tan(/3) and 
tan(0/) terms, which in (6.5) are not included. We refer to [28] for a review of possible lens 
equations. 
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Figure 12. (a): for the Virbadhra-Ellis lens equation in the neutrino case, for a black hole 

with M = 10® Mq and with Dol=^^ Kpc, Z9os=19 Kpc. The numerical solution for the classical 
and the energy-dependent result, (b): /3(0/) as in (a) but for a 1 GeV neutrino beam. 


7 Radiative effects and the geometry of leasing 

Turning to our case study, radiative effects in the lens equations can be introduced by 
replacing the expression of the angular deflection generated by the source on the source 
plane, which is a function of the impact parameter b (a = a{b)) with the new, energy 
dependent relation a{b, E) whose general form is given by (4.9). 

For simplicity we consider a pointlike source, and a pointlike deflector, as shown in Fig. 
10. We recall that for a massless particle the geodesic motion is determined in terms of 
the energy E and of the angular momentum L at the starting point of the trajectory. The 
gravitational deflection, however, can be written only as a function of the impact parameter 
b of the source, with b = E/L, which is an important result of the classical approach. For 
a further clarification of this aspect, which differs from the semiclassical analysis we are 
interested in, we briefly overview the classical case, using the lens geometry as a reference 
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point for our discussion. 

For a source located on the source plane at an angular opening j3 (in the absence of the 
deflector), the initial conditions can be expressed in terms of the two components of the 
initial momentum p = (pr,P(f>) on the plane of the geodesic, or, equivalently, by the pairs 
{pr,E) or {p^, E), with E the initial energy of the beam. We recall that for a Schwarzschild 
metric these are defined as 

2 ^ 2GM\ . 2 ■ 2 

r, p^ = -r Pt = 1 1-I P<p = -r sm ifcp. 

(7.1) 

We have denoted with x = dx/ds the derivative respect to the affine parameter, pt and p^ 
related to the energy and to the angular momentum as pt = E and p^ = —L, and with the 
motion taking place on the plane i? = 7r/2 (p^ = 0). They are constrained by the mass-shell 
condition 




r J 



T J 


-r'^{p^f = 0 , 


with {p^ = r,p^ = (^,p^ = i). 

The lens equation, usually written as 


(7.2) 


L{/3,9i)=0, (7.3) 

can also be written, equivalently, in the form of a constraint between fd and h using (6.19). 
We can use any of the independent variables mentioned above. For a given initial mo¬ 
mentum of the beam, emitted from the plane of the source, the lens equation will then 
determine the position of the source in such a way that the geodesic motion will reach the 
observer at its location on the optical axis. In particular, an interesting description emerges 
if we choose as initial conditions the angular position of the source {fd) and the value of the 
impact parameter b. These two conditions fix the direction of the trajectory of the beam at 
its origin on the source plane. In these last variables, the lens equation will then determine 
one of the two in terms of the other in such a way that outgoing geodesic will reach the 
observer. 

The inclusion of an energy dependence in the angle of deflection a renders this picture 
slightly more complex. For instance, the lens equation will now depend on 3 parameters, 
which can be chosen to to be {/d,6i,E) or {fd,Pr,P(f)) or any other equivalent combination, 
with one of the three fixed in terms of the other two by the equation itself. For a monochro¬ 
matic and spherical source of energy E, fixed at a position /?, emitting a beam with a given 
impact parameter b respect to the deffector, the lens equation may not have a real solution, 
since the deflector may disperse the beam in such a way that it will never reach the observer. 
For a fixed spherical source which emits photons or neutrinos of any energy, one can look 
for solution in the reduced variables b, E. Being b related to the primary and secondary 
images 9i±, the beam that reaches the observer will be characterized by a unique energy 
E^ assuming that the images are detected at angular positions 9j±. 

The argument above can be repeated by using any triple combination of independent kine¬ 
matic variables among those mentioned above. 
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Having clarified this point, we now move to a description of the actual implementation of 
the lens equation is this extended framework. The angular location of the image Oj and the 
impact parameter are related in the geometry of the lens by Eq. (6.19), and this allows to 
search for solutions of the lens equation (6.18) in regions characterized by smaller values of 
the impact parameter (20 < < 100) where the angular deflections are stronger. 

The key to the derivation of the radiative lens equation are Eqs. (4.9) and (6.19). Com¬ 
bining the two relations we obtain 
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where the ellipsis refer to the extra logarithmic contributions present in Eq.(4.9). The 
expression above is known analytically if we manage to solve explicitly the semiclassical 
equation (4.1), otherwise it has to be found by a numerical fit. However, it is clear that the 
ansatz for the fit has, in any case, to coincide with Eqs. (4.9) and (7.4), due to the typical 
functional forms of the solutions of Eq. (4.1). For instance, in the case of a thin lens, the 
modifications embodied in (7.4) can be incorporated into the new equation 


/3 = 9j-a{b{9i,E))-:^, (7.5) 

Dos 

which is an obvious geueralizatiou of (6.9), the latter being valid ouly in the classical GR 
case. As we are going to illustrate below, (7.5) can be studied numerically for several 
geometrical configurations, which are obtained by varying the lensing parameters D^s und 
Dol- 

A similar approach can be followed for the VE or for any other classical lens equation. The 
insertion of a{9i,E) given by (4.9) into (6.18) generates the radiative lens equation 

Dos tan/3 = Dos tan 9j — + tan(a(0z, E) — 9j)), (7-6) 

which takes into account also the quantum corrections and is now, on the contrary of (6.18), 
energy dependent. At this point it is clear that all the lens observables, such as magui- 
ficatious, shears, light curves of microlensiug etc. desceud rather directly by this general 
prescription. 

For instance, we can determine for the Virbadhra-Ellis lens the expression for the magnifi¬ 
cation using the radiative (semiclassical) expression 

/a = — (7.7) 

X2 

X 1 = Dossin9i (1 -I- ((Hoitan^/ -|- {Dol - Dos)tan{a{9i, E) - 9i))/ Dos)‘^)^^‘^, 

X 2 = (T>oLtan0/ -I- {Dol - Hqs) tan(a(0/, i3) - 0/)) 

x(sec^ 9i -I- {Dol — Dos)/Dos{sec^ 9i + sec^{a{9i, E) — 9i) {a'{9j, E) — 1))), 
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Figure 13. /3(0/) with 20 < bh < 100 in the neutrino case for a black hole with M = 4.31 10® M©: 
D^s = 1 kpc (a), Dls = 1 pc (b), D^s = 1 milliparsec (mpc) (c). 


where a' = dajdOi. As clear from Eqs. (7.6) and (7.7), both equations are very involved, 
although they can be investigated very accurately at numerical level. It is also possible to 
discuss the analytical form of the solutions within the formalism of the 1 /6” expansion. In 
fact, we are entitled to expand all the observables of the fully nonlinear lens in the angular 
deflection a, and work at a certain level of accuracy in the angular parameters. In this 
work, however, we prefer to proceed with a direct numerical analysis of the full equations, 
both for the thin and for the VE lens, leaving the discussion of the explicit solutions to a 
future work. 

7.1 Numerical analysis for neutrino and photon lensing 

The analysis that we present in this section, especially in the neutrino case, is of exploratory 
nature. It has the aim to test the consistency of the theoretical approach presented in the 
previous sections, rather than being an explicit proposal for the detection of such effects. 
We start investigating the behaviour of the solutions for the VE lens equation, which are 
shown in Figs. 12, by plotting the angular position of the source (/3) as a function of the 
location of the image 6i, for neutrino beams. In panel (a) we have chosen distances between 
lens/source and observer of galactic size. The branch of the solution with Oj > 0 describes 
a primary image, while for Oj < 0 the points (/3, Oj) on the curve describe configurations 
corresponding to secondary images. One can test the energy dependence of the solution by 


29 - 













Dls = I Kpc 


Dls = 1 pc 


-1.x 10-'^ 
-2.xl0‘‘^ 
-3.xl0‘‘^ 
-4.xl0"‘^ 
-5.xl0"‘^ 
-6.xl0"‘^ 
-7.xl0‘‘^ 


0.2 0.4 0.4, 0.8 1.0 


Classical 
E = 1 GeV 
E = 1 TeV 
E = 1 PeV 


(a) 


-l.xlO-® 

-2.X10-® 

-3.X10-® 

-4.X10-® 

-5.X10-® 

-6.X10-® 


0.2 0.4 0,6 0.8 1.0 


(b) 


Ols = 1 mpc 



(c) 


Figure 14. (Sagittarius) with 20 < bh < 100 in the neutrino case for a black hole with 

M = 4.31 X 10® M©: D^s = 1 Kpc (a), D^s = 1 pc (b), Dls = 1 nipc (c). The strong suppression 
of the magnification parameter for panels (a) and (b) are associated with the secondary image. 


varying the energy of the original beam in the lens equation (7.6), which, for this specific 
geometry, is clearly unnoticeable. In general, in fact, specific geometries {Dql, Dls) select 
impact parameters bh in the beam path which are quite large. In this case the impact 
parameter turns out to be pretty large {bh ~ 10^), corresponding to very weak deflections, 
and causes a superposition between the two curves, the one describing the classical GR 
solution of (6.18), and the semiclassical one, obtained by solving (7.6). The curves intersect 
the 9j (image) axis in two opposite points {9i = ±0^;), giving rise to the Einstein ring, 
which are obtained for /3 = 0, i.e. for a complete alignment of the lens/source/observer 
along the optical axis. In this figure the negative {3 range is symmetric and hence it is not 
shown. It can be obtained by a parity flip of the two positive branches with [3 —>■ —j3 and 
91 —> -9i. 

In Fig. 12 (b) we investigate the dependence of the classical solution for three values of Dls- 
The solution curves so generated exhibit variations which are clearly far more significant 
that any radiative correction which might affect the lensing geometry. Notice that for a 
given solution of the lens equation 9j, as we increase the distance Dls, for a given angular 
position of the image, the source moves towards the optical axis. This behaviour, in the 
secondary image, appears to be reversed: in this case larger values of Dls require larger 
angular values of /3, for a given angular position 0/ of the image. 
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Figure 15. /3(0/) with 20 < bh < 100 in the neutrino case for a black hole with M = 1.40 10® M©: 
= 1 Kpc (a), Dls = 1 pc (b), Dls = 1 nipc (c) (Andromeda galactic center). 


Solutions with smaller values of the impact parameters are those which are more favourable 
from the point of view of the semiclassical treatment, since in these cases the deflections 
are larger and induce larger gradients into the lens equation (7.6). For this reason we 
investigate two leasing configurations corresponding to the case of the supermassive black 
holes located at the center of our galaxy (Sagittarius A*, 8 kpc) and of the nearby galaxy 
Andromeda {Dql ~ 780 kpc), and vary the distance Dls between the source and the lens. 
We show in Fig. 13 the solutions for the Sagittarius configuration. In panels (a) and (b) the 
solutions that we identify correspond to secondary images obtained for very small values 
of 0/, of the order of a milliarcsecond. These are the only configurations which guarantee 
close encounters between the cosmic ray beam and the black hole, with 20 < bh < 100. As 
the distance Dls gets reduced, the equation has solutions with values of Oj which define 
primary images, being /? > 0 (panel c). At the same time, as shown in panel (c), the angular 
position of the source moves from ,0<Oto/3>O. A primary image is shown to form when 
Dls is 1 milliparsec (mpc). 

Simple considerations show that such lensing configurations are not unreasonable. For 
instance, for the supermassive black hole that we consider (with M=4.31 x 10®Mq), this 
distance is of the order of 5 x lO^rs, with denoting its Schwarzschild radius (~ 6 x 10® 
km). Being the center of our galaxy rather densely populated by massive compact sources. 
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Figure 16. ^{Oj) with 20 < bh < 100 in the neutrino case for a black hole with M = 1.40 10® Mq: 
D^g = 1 Kpc (a), Dgg = 1 pc (b), Dgg = 1 mpc (c) (Andromeda galactic center). The strong 
suppression of the magnification parameter, in this case, is associated to the secondary image for 
this geometric configuration. 


one could envisage a distribution of these covering a large array of possible distances from 
the center of the black hole. For instance, it has been found that stars may orbit the 
supermassive galactic black hole with orbital periods even of the order of T ~ 11.4 ys, 
corresponding to orbital distances (R) from its center as close as few astronomical units 
(AU), R = ~ 5 AU, i.e. R ~ 125 r^. While such distances may correspond to realistic 

lensing configurations, the lensing resolutions of these specific events, which is of the order 
of a milliarcsecond, remains a challenging aspect of these studies. This is due to the strong 
quasi-alignment required between source, lens and observer along the optical axis, which 
might be difficult to measure. 

The energy dependence of the lensing configuration, extracted from Eq. (7.6), is illus¬ 
trated by plotting the solutions for several values of the initial energy of the neutrino beam, 
corresponding to 1 GeV, 1 TeV and 1 PeV respectively. A comparison with the classical 
GR solution is included. Differences among the various predictions for the position of the 
source can be as large as 15%. 

The analysis for the magnification of this lensing configuration is presented in Fig. 14. Pan¬ 
els (a) and (b) show the strong suppression of the secondary images identified in Fig. 13 (a) 
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Figure 17. /3(6*/) with 20 < bh < 100 in the photon case for a black hole with M = 4.31 10® M©: 
Dls = 1 Kpc (a), Dls = 1 pc (b), D^s = 1 mpc (c) (Sagittarius). 


and (b). The magnification of the secondary and primary images of Fig. 13 (c) is shown 
in Fig. 14 (c). The /r < 0 and /r > 0 regions, in this figure, are separated by asymptotes 
for 9 1 = Be, the Einstein angle of the lens. These regions in /r correspond to the secondary 
(/3 < 0) and (/3 > 0) branches of the panel (c) in Fig. 14 and therefore refer to secondary 
and primary pictures respectively. 

The dependence on the energy of the incoming neutrino beam appears in the form of 3 
displaced curves and respective asymptotes. The three vertical asymptotes therefore char¬ 
acterize the dependence of the Einstein radius on the energy. 

The analysis is repeated in Figs. 15 and 16 for lensing events detected on Earth from the 
galactic center of the Andromeda galaxy. In this case as before, we vary the distance be¬ 
tween the lens and the source, with Dls = ^ kpC) 1 pc and 1 milliparsec in panels (a), (b) 
and (c) respectively. The first two plots correspond to secondary images while panel (c) 
describes a primary image solution of the lens equation. The corresponding plots for the 
magnifications, in the three cases, are shown in Fig. 16, evidencing its strong suppression 
for the two secondary images (panels (a) and (b)), and the enhancement for the primary 
image in Fig. 15 (c), shown in panel (c). Also in this case we illustrate the dependence of 
the image solutions on the energy of the incoming neutrino beam. 

The patterns found for the neutrino lenses remain valid also for the photons, as one can 
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Figure 18. /i(0/) Sagittarius with 20 < bh < 100 in the photon case for a black hole with M = 4.31 
10 ®Mq: Dls = 1 Kpc (a), D^s = 1 pc (b)) -Dls = 1 mpc (c) (Sagittarius). 


easily figure out by a cursory look at Figs. 17 and 18 for the supermassive black hole in 
Sagittarius. In panels (a) and (b) of Fig. 17 the secondary images found as solution of 
the lens equation are suppressed in magnitude, as shown in Fig. 18 (a) and (b), while the 
solution in Fig. 17 (c), corresponding to one primary and one secondary image, is associ¬ 
ated with the magnihcation given in Fig. 18 (c). In these two sets of hgures the energy 
dependence of the result is quite small. 

Finally, the numerical result for the photon lensing case, generated by the supermassive 
black hole in the Andromeda galactic center is discussed in Figs. 19 and 20. While the 
two secondary images found in 19 (a) and (b) give suppressed magnifications, the solution 
in panel (c) corresponds to a primary image. The corresponding magnihcation, shown in 
Fig. 20 has, obviously, a single branch, with an Einstein radius located at approximately 
7 x 10“® arcsec. 


8 Post-Newtonian corrections: the case of primordial black holes 

We have seen in the previous sections that the bh(ct) expression for the dehection does not 
suffer from any apparent divergence (from the gravity or external held side) due to well- 
dehned structure of the Newtonian cross section. The expression given in (3.12), in fact, is 
similar to the ordinary Rutherford scattering encountered in electrodynamics. 
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Figure 19. /3(6*/) with 20 < bh < 100 in the photon case for a black hole with M = 1.40 10® M©: 
Uls = 1 Kpc (a), Dls = 1 pc (b), Dls = 1 mpc (c) (Andromeda). 


The dependence of the resulting cross section on the scale GM/cP, the Schwarzschild radius, 
manifests as an overall dimensionful constant. Therefore, the inclusion of the electroweak 
corrections - and the logarithmic dependence on the energy of the terms in the expansion 
that follows - do not appear in combination with the macroscopic scale This allows, in 
principle, an extension of the perturbative computation up to any order in the electroweak 
coupling constant a^. It is also clear that this result is expected to be valid for any 
renormalizable held theoretical model, when combined with an external static gravitational 
held of Coulomb type, as in the case of the Newtonian limit of GR. 

From now on, we will be using the notation nPN to indicate the (post-Newtonian) order in 
the potential at which we expand the Schwarzschild metric. For instance, contributions of a 
certain nPN order involve corrections in the external held proportional to <1)”+^, with OPN 
denoting the ordinary (lowest order) Newtonian (i.e. zeroth post-Newtonian) contributions 
proportional to <h, as given in Eq. (3.7). The inclusion of the higher order corrections 
in the external potential modihes this simple picture due 1) to the need of introducing a 
cutoff regulator in the computation of the Fourier transform of the higher powers of the 
Newtonian potential and 2) to the presence of the Schwarzschild radius in the actual 
expansion. These features emerge already at the hrst post-Newtonian order (IPN) for an 
uncharged black hole and at order OPN for the Reissner-Nordstrom (RN) metric (charged 
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Figure 20. /i(0/) with 20 < bh < 100 in the photon case for a black hole with M = 1.40 10® M©: 
Dls = 1 Kpc (a), Dls = 1 pc (b), Dls = 10 mpc (c). The strong suppression of the magnifica¬ 
tion parameter, in this case, is associated to the secondary image for this geometric configuration 
(Andromeda). 


black hole). 

Both points 1) and 2) are, in a way, expected, since the microscopic expression for the 
transition matrix element given by (3.1), in fact, cannot be extrapolated to the case of a 
macroscopic sonrce, with the presence of a macroscopic scale snch as the black hole horizon. 
This seems to indicate that the success of the Newtonian approximation is essentially dne 
to the rescaling of fonnd in the expression of the cross section, which is a featnre of 
this specific order, and is therefore limited to a 1/r potential. It is then natnral to ask if 
there is any other realistic case in which the post-Newtonian corrections can be inclnded 
in an analysis of this type. Obvionsly, the answer is affirmative, as far as we reqnire that 
Vs is microscopic and that the energy of the beam, which is an independent variable of 
a scattering event, is at most of the order of l/r*. Under these conditions, we are then 
allowed to extend onr analysis throngh higher orders in with scatterings in which the 
dimensionless parameter VgQ with q the impact parameter, is at most of 0{1). This specific 
sitnation is enconntered in the case of primordial black holes, where Tg can be microscopic. 
We are going to illnstrate this point in some detail, since it becomes relevant in the case of 
primordial black holes. 
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8.1 Post Newtonian contributions in classical GR 


To illustrate this point we extend the expansion of the Schwarzschild metric at order OPN 
given in (2.2). A similar expansion will be performed on the RN metric. 

For this purpose, it is convenient to perform a change of coordinates on the Schwarzschild 
metric 


ds^ 



2GM\ 
r J 


dt^ 



2GM 


-1 


dG 


r^dn 


( 8 . 1 ) 


in such a way that this takes an isotropic form. The radial change of coordinates is given 
by 


which allows to rewrite (8.1) as 




( 8 . 2 ) 


ds'^ = A{p)dt^ — B{p){dp^ + p^ dfl) 


(8.3) 


with 


A{p) = 


(1 - GM/2pf 


B{p) = {l + GM/2pY 


(8.4) 


{l + GM/2pf 

Post-Newtonian (weak field) corrections can be obtained by an expansion of A and B taking 
M/p <C 1. Up to third order in this is given by 


A{p) = 1 + 2^ + 2^"^ + 

B{p) = 1 - 24> -h ^$2 - 


(8.5) 

( 8 . 6 ) 


In the RN spacetime for a charged black hole the analysis runs similar. The interest 
in this metric is due to the fact that the lowest order potential, in this case, involves 
charge-dependent 1/r^ contributions which, for an uncharged black hole, appear at first 
post-Newtonian order (IPN). The metric, in this case, is given by the expression 
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ds^ =\l- 
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(8.7) 


with Q denoting the overall charge of the black hole. It has two concentric horizons which 
become degenerate in the maximally charged case. The two horizons are the solution of the 
equation 


2GM _ 

T o 1 0 

T } 


( 8 . 8 ) 


with solutions r = GM ± — GQ"^. The RN black hole has a maximum allowed 

charge Q = M\/G, in order to avoid a naked singularity. In this case, the radial change of 
variables which brings the metric into a symmetric form is given by 


r = p 


^ gm + Vgq ^ 


^ gm-Vgq ^ 


(8.9) 
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so that the RN spacetime in isotropic coordinates is 
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1 - 
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( 8 . 10 ) 


We just recall that for a massless particle in this metric background the angle of deflection 
and the impact parameter are given by the expressions 
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( 8 . 12 ) 


where tq is the closest distance of approach. It’s convenient to normalize r, ro and Q to 
the Schwarzshild radius r, = 2 GM and introduce the variables 


X = 


Xq = 


ro 


d = 


Q 


2 GM 2 GM " 2 GM 

With this redefinitions the deflection can be expressed in the form [29] 

a(a;o) = G(xo) F((/)o, A) - vr 


with 


G{xo) = 


4x0 


./l - ^ VVi - rV){r2 - rV 

V ^0 Xq 

and with 

rpo 

F((/)o,A)= / {1 — X sitV (j))~^^'^d(j) 

Jo 

being an elliptic integral of the first kind with arguments 


(8.13) 
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(po = arcsin 
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The ri are the roots of the fourth order polynomial 
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Figure 21. Comparison of the deflection angle for the Schwarzshild case and the maximally charged 
Reissner-Nordstrom case in the near-horizon (a), in the very-near-horizon (b). 


ordered so that ri > r 2 > > r 4 ^. The comparison between Schwarzschild and RN 

deflection angle is shown in Figure 21. The plots describe the behaviour of the angular 
deflection as a function of the impact parameter for a RN and Schwarzschild metric in 
the region with 10 < 6^ < 50 (top left) and 4 < bh < W (top right) for the maximally 
charged case. The differences tend to be very pronounced as we approach the horizon of 
the Schwarzschild metric. 

As pointed out in [24] in the Schwarzschild case, the 1/6 expansion for the deflection angle 
does not reproduce the photon sphere singularity of the Schwarzschild metric, which is 
achieved using the exact GR expression in terms of elliptic function given in (8.14), but 
it represents nevertheless an improvement respect to the OPN order. Expanding the RN 
metric in M/p 1 up to the third order, the 2PN approximation gives 
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( 8 . 21 ) 


Inserting this expansion into the deflection integral, we can account in a systematic way of 
the 1 /h corrections in the angle of deflection a 


GM / GQ^\37rG^M^ 

a(l ,)=4 —+ 5-3^ + 
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( 8 . 22 ) 


The deflection (8.22) in the maximally charged case is given by the expression 


„, GM ^ {GM f 80 (GM)3 
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62 
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(8.23) 


In the next subsection we are going to illustrate how the inclusion of these expansions at nPN 
order affects the computation of the quantum corrections to the angular deflection. The 
corrections are embodied in a geometric form factor whose expression is entirely controlled 
by the 1/6 expansion. 
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8.2 Quantum effects at 2nd PN order 


The inclusion of the PN corrections to the external background requires a recalculation of 
the cross section, with the inclusion of the additional terms in the fluctuation of the metric 
in momentum space. As usual we consider a static source, so that the metric is written as 


hf^uiq) = 2TT5{qo)hf,„{q). 


( 8 . 24 ) 


At leading order in the external held <1> both the timelike and the spacelike components are 
equal ( /iqo = ha), while at higher orders they are expressed in terms of two form factors 
ho and hi 

h^u{q) = ho{q)5ofiSou + hi{q){rj^^ - ( 8 . 25 ) 

which at higher order in the weak external held are given by 
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where we have explicitly reinstated the dependence on the speed of light. Below we will 
conform to our previous notations in natural units, with c = 1. 


• Neutrinos 


The computation, at this stage, follows rather closely the approach of the previous sections, 
giving for the averaged squared matrix element in the neutrino case 




i>2 (Pl > P2 )|^ 1 (pi, P2 ) 


{q)hpAQ), ( 8 . 27 ) 


with T being the time of the transition, and the differential cross section 
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dW 
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duf. 


( 8 . 28 ) 


We have denoted with duf the density of hnal states in the transition amplitude, and with 
ji the incoming hux density. After integration over the hnal states, and using \pi\ = \p 2 \, 
we obtain the expression 
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( 8 . 29 ) 


where we have introduced the gravitational form factor of the external source 

Fg{q) = (ho{q) - hi{q)J . ( 8 . 30 ) 

Notice the complete analogy between the corrections coming from a distributed source 
charge, for a potential scattering in quantum mechanics, and the gravity case. In the 
evaluation of Fg in momentum space we are forced to introduce a cutoff A, being the 
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Fourier transforms of the cubic contributions in divergent. The singularity is generated 
by the integration around the region of r ~ 0 in the Fourier transform of the potential. The 
relevant integrals in this case are given by 

1 


with 


In = (VX 


dvr 

-^1 -*2 ’ 
qz 


iq-x 

(8.31) 

27r2 

(8.32) 



and with Is requiring a regularization with an ultraviolet cutoff in space (A) 

sin(|^|r) 


^ dvr 

l^l , 


dr 


(8.33) 


The choice of A is dictated by simple physical considerations. Given the fact that consis¬ 
tency of the expansion requires that h is clear the appropriate choice in the 

regulator is given by the condition that this coincides with the Scwarzschild radius, i.e. 
A ~ r^. Expressed in terms of the cutoff, we obtain for the geometric form factors the 
expressions 
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where we have indicated with Ci the cosine integral function 

Ciw= 

J oo ^ 

From the previous equations we obtain the cross section 


-|g|Ci(A|g|) (GM)= 
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-|g|Ci(A|g|)J(GM)3j , 

(8.36) 

which is valid at Born level and includes the weak field corrections up to the third order 
in <h. In the expression of the cross sections, we use the subscript nPN, with n = 0,1,2 
to indicate a n-th order expansion of the metric in the gravitational potential, while the 
superscripts ((0), (1) and so on) label the perturbative order in a^. The leading order cross 
section at order 2PN, for instance, takes the form 


da 

dQ 


( 0 ) 


2PN 
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rM2{E,e), 
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with 


VM2{E,9) = 


1 - ^{GM)Es\n^- + ^(GM)2 Esin sin 


2 A E sin - 
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Q 

— 2E sin - Ci ( 2 A E sin - 


(8.38) 


dl 



where we have factorized the tree level result du/dniQpj^ given in (3.12). The post- 
Newtonian form factor VJ\f 2 {E, 9) induces an energy dependence of the cross section which 
is unrelated to the electroweak corrections. The analysis, in fact, can be extended at one 
loop in the electroweak theory. In this case, a lengthy computation gives the 2PN result 


da 


( 1 ) 


dVt 


2PN 


da 

dQ 
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AGf 




VAf2{E,e), 

(8.39) 


where we have inserted the one loop expression given in (3.19). 

We can obtain an explicit solution of the corresponding semiclassical equation at order IPN. 
Using the expression of the VAf function at this order 


VA/i{E,e) 



TT 9 

-{GM)Esin- 
8 2 


(8.40) 


on the right hand side of (8.38) in order to generate the IPN cross section at Born level, 
and solving the corresponding semiclassical equation (4.1) we obtain 


^^liPN ~ ^ {GM)‘^ ^ — 1 -|- csc^ ^ 21nsin + E (GM)^7r^4 -|- (cos a — 3) esc 

- E‘^{GM)'^'k^(\ + cos a + 4 In sin —'j. 

32 V 2 / 


(8.41) 


At this point, we can invert Eq. (8.41) for a{b) obtaining 
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( 0 ) _ 




hH 
37r^ 


-^E\GMf)+0{bi) 


(8.42) 


for the tree level post Newtonian one. 

For the Reissner-Nordstrom geometry the situation is similar. The post-Newtonian form 
factor is then given by 


VM{E, ( 


\RN 
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(8.43) 


and the impact parameter in the IPN approximation is 

= 4 {GM)‘^(^ -1 + csc^ ^ -h 21nsin + E {GM f[l + 3 ^^2 )^ + (cos^ - 3) esc 
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The inversion formula in this case is 
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(8.45) 


• Photons 

We can extend the analysis presented above for neutrinos to the photon case. Here the 
cross section takes the form 
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and, as in the neutrino case, we have 
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where we inserted the tree level cross section for the photon 
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(8.48) 


In the OPN Newtonian limit, this cross section has been computed in [5], and takes the 
form 
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and 
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are the relevant electroweak form factors entering in the computation. In the previous 
equations the sum fk is over all Standard Model fermions, with rrifk and Qfj^ their masses 
and charges. is 1 for leptons and 3 for quarks. Proceeding similarly to the neutrino 
case, the one loop cross section in the 2PN approximation takes the form 


da 
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7,2PN 


da 


( 1 ) 

FAA2(F,0), 

7,0PN 


(8.52) 


with VN '2 given by (8.38), which can be inserted again in (4.1) and investigated numerically. 
Solving at order IPN the analogous of (8.52), the solution of (4.1) gives 
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V ^ ^/o V ^/V4 4/ 


- E^(GM)^7r^ (11 + 12 cos a + cos 2a + 32 In sin —) . 

256 V 2/ 

In the photon case the inversion formulae at orders OPN and IPN are given by 


(8.53) 


“l^loPN ^ ^ ~ ^ ~ 1) ~ ^ 


and 
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li'ipN - 4 

-li£(GAC-l)+0(6j) 


(8.55) 


respectively. 


8.3 Range of applicability 

The structure of the one-loop 2PN result for neutrinos and photons shows the complete fac¬ 
torization between the quantum corrections and the background-dependent contributions. 
While the former are process dependent, the latter are general. Obviously, this result is not 
unexpected, and follows rather closely other typical similar cases in potential scattering in 
quantum mechanics. An example is the case of an electron scattering off a finite charge 
distribution characterized by a geometrical size R, where the finite size corrections are all 
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contained in a geometric form factor. 

We recall that for a Coloumb interaction of the form V (r) = ^ jr, the cross section is given 
in terms of the pointlike (p) amplitude 



,2 


(8.56) 


with q = k — k' and = 2\k\ sin0/2 being the momentum transfer of the initial (final) 
momentum of the electron k {k') and charge e. The scattering angle is measured with 
respect to the z-direction of the incoming electron. The charge of the static source has also 
been normalized to e. The corresponding cross section is given by 


=\f(Q\ |2 _ (2m)^e^ 

drip ^ 16/c^sin^0/2 


(8.57) 


and the modification induced by the size of the charge distribution {p{x)) is contained in 


m 


/ 


dxp{x)e^’^'^ 


(8.58) 


with 



(8.59) 


For a uniform charge density, for instance, the geometrical form factor F(q), which is the 
transform of the charge distribution, introduces a dimensionless variable qR in the cross 
section which is absent in the point-like (Coulomb) case, of the form 



(8.60) 


The validity of the expression above is for qR < 1, and the presence of the geometrical 
form factor is responsible for the fluctuations measured in the cross section as a result of 
the finite extension of the charged region. 

In the analysis of the nPN corrections in gravity, the situation is clearly analogous, with 


the size of the horizon taking the role of the classical charge radius R. For ordinary (macro¬ 


scopic) horizons (e.g. of a km size) ~ GM invalidates the perturbative expansion due 
to the appearance of the GME parameter in the expression of the post Newtonian factor 
0), which is small only if FI ~ 1/GM, a choice which is not relevant for our analysis, 
since it applies to particle beams whose energy is in the very far infrared. 

By imposing that the cutoff A coincides with the Schwarzschild radius (A = r®), one 
can immediately realize that the post-Newtonian expansion gets organized only in terms 
of this parameter (GME). In the regions of strong deflections, which are those that con¬ 
cern our analysis, we can reasonably assume that y = sin0/2 ~ 0(1), if we use the GR 
prediction to estimate the bending angle. This allows to discuss the convergence of the 
PN expansion only in terms of the energy E of the incoming beam and of the size of the 
horizon. The analogous of the charge oscillations given by (8.60), in the gravitational case, 
are then uniquely related to the post-Newtonian function VM, and hence to the size of the 
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Figure 22. Comparison of nPN approximations for a{b) in the photon case with Mpbh = 
10“^®Mq and for EGM « 1 (a). In (b) and (c) we show the VAf function for different ener¬ 
gies. 


parameter AE ~ VgE which defines its expansion in powers of the gravitational potential. 
Assuming a small value of x = GME, we can indeed rewrite (8.38) via a small-x expansion, 
obtaining 

VM{x, y) = 1- ^xy + x^y"^ (1 - 7 ^; - log x - log 4y). (8.61) 

O 

This expression can be used to investigate the range of applicability of these corrections in 
terms of the two factors appearing in x, the energy of the incoming beam and the size of 
the horizon of the gravitational source. The requirement that such a parameter be small 
defines a unique range of applicability of such corrections in the quantum case. 

One possible application of the formalism which renders the PN corrections to the gravi¬ 
tational scattering quite sizable is in the context of primordial black holes [30], which have 
found a renewed interest in the current literature [31, 32]. 

We just mention that primordial black holes (PBHs) have been considered a candidate com¬ 
ponent of dark matter since the 70’s, and conjectured to have formed in the early universe 
by the gravitational collapse of large density fluctuations, with their abundances and sizes 
tightly constrained by various theoretical arguments. These range from Hawking radiation, 
which causes their decay to occur at a faster rate compared to a macroscopic black hole (of 
solar mass); bounds from their expected microlensing events; their influence on the CMB, 
just to mention a few [33]. For instance, the mechanism of thermal emission by Hawking 
radiation sets a significant lower bound on their mass (~ 5 x 10 ^^( 7 ), in order for them to 
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survive up to the present age of the universe. This bound satisfies also other constraints, 
such as those coming from the possible interference of their decay with the formation of 
light elements at the nucleosynthesis time. With the launch of the FERMI gamma ray space 
telescope [11], the interest in this kind of component has found new widespread interest. 
The unprecedented sensitivity of its detector in the measurement of interferometric patterns 
generated by high energy cosmic rays (femtolensing events), such as Gamma Ray Bursts 
[10], has allowed to consider new bounds on their abundances [34]. The hypothesis of hav¬ 
ing PBHs as a dominant component of the dark matter of the universe provides remarkable 
constrains on their allowed mass values, except for a mass range 10^®kg < Mpbh < lO^^kg, 
where it has been argued that they can still account for the majority of it. In other mass 
ranges several analyses indicate that the PBH fraction of dark matter cannot exceed 1% of 
the total [33]. 

PN corrections turn out to be significant for PBH in this mass range, due to the large 
variation induced on the VM function by the IPN and 2PN terms. These may play a 
considerable role in a PBH mediated leasing event. We illustrate this behaviour by showing 
plots of the post Newtonian behaviour of the relevant expressions for leasing. In Fig. 22 
(a) we plot the angular deflection as a function of the impact parameter for the Newto¬ 
nian OPN, and relative post Newtonian corrections. We have considered a primordial black 
hole with a mass of 10“^® Mq-, which carries a microscopic Schwarzschild radius (300 fm) 
and chosen E = 1/{GM) = 0.6 MeV for the incoming photon beam. The impact of the 
corrections on the gravitational cross section are quite large, as one can easily figure out 
from panel (b), where we plot the factor VAf as a function of the Schwarzschild radius for 
these compact massive objects, for ~ 1 fm. For a more massive primordial black hole, 
with 200 < bh < 1000, the pattern is quite similar, as shown in panel (c). In both cases 
the post-Newtonian corrections appear to be significant, of the order of 15-20 % and could 
be included in a more accurate analysis of leasing for these types of dark matter candidate 
solutions. 


9 Conclusions and perspectives 

We have presented a discussion of neutrino leasing at 1-loop in the electroweak theory. In 
our approach the gravitational field is a static background, and the propagating matter 
fields are obtained by embedding the Standard Model Lagrangian on a curved spacetime, 
as discussed in previous works [14, 17]. As in a previous study [5], also in our current case 
the field theoretical corrections to the gravitational deflection are in close agreement with 
the predictions of general relativity. The agreement holds both asymptotically, for very 
large distances from the center of the black hole, of the order of 10® horizon sizes (6/i), but 
also quite close to the photon sphere (~ 20 bh)- In this respect, the similarity of the results 
for photons and neutrinos indicates the consistency of the semiclassical approach that we 
have implemented. As noticed in [19], the inclusion of the quantum effects causes the 
appearance of an energy dependent dispersion of a particle beam, which implies a violation 
of the classical equivalence principle. 
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Various types of lens equations have been formulated in the past using classical GR, 
and we have illustrated the modifications induced on their expressions by the inclusion of 
the suppressed l/fe” corrections in the impact parameter to the angular deflections. We 
have then developed a formalism which allows to include the semiclassical results, due to 
the radiative effects in the propagation of a photon or a neutrino, in a typical lensing 
event. We have considered both the case of a thin lens, which is quadratic in the deflection 
angle, and the fully nonlinear case, taking as an example the Virbhadra-Ellis lens equation. 
Radiative and post-Newtonian effects induce a dependence of the angle of deflection with 
the appearance of extra l/ft” suppressed contributions and of extra logarithms of the impact 
parameter, that we have studied numerically for some realistic geometric configurations. In 
general, radiative effects are significant only for configurations of the source/lens/observer 
which involve small impact parameters in the deflection (bh ~ 20), and require angular 
resolutions in the region of few milliarcseconds. Our results are valid for a Schwarzschild 
metric, considered both in the Newtonian and in the post-Newtonian approximation, but 
they can be extended to other metrics as well. 

We have also discussed the consistency of the post-Newtonian approach. We have shown 
that such corrections can be consistently taken into account in the case of microscopic 
horizon sizes, such as primordial black holes. These corrections have been shown to factorize 
and be accounted for by a post-Newtonian function. Our analysis can be extended in several 
directions, from the case of Kerr-Newman metrics to the study of microlensing and Shapiro 
delays, and to dynamical gravity. We hope to return on some of these topics in a future 
work. 
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A Fluctuations to first order in the Newtonian potential 

The solution for the static massive source is obtained from the linearized equation 

□ (v - (A.l) 

where h = h^’^r]^iy, and can be rewritten as 

ah^u = KS^„ with 5^, = - . (A.2) 

The external field hf^u is obtained by convoluting the static source with the retarded prop¬ 
agator 

n ( \ 1 *^(^0 -yo-\x-y\) , , 

GR{x,y) = - --—--, (A.3) 

dvr \x — y\ 
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normalized as 


(A.4) 


n\GR{x,y) = S^{x - y). 


The solution of Eq.(A.2) takes the form 



(fy GR{x,y) S^y{y) 


(A.5) 


with the EMT of the external localized source, dehning given by 



(A.6) 


For a compact source of mass M at rest at the origin, with = (M, 0), we have 



(A.7) 





(A.8) 




(A.9) 


where the held generated by a local (point-like, L) mass distribution has the typical 1/r 
(r = |x|) behaviour. The huctuations are normalized in such a way that hfj^y has mass 
dimension 1, as an ordinary bosonic held, with n of mass dimension —1. 

B l/b^ corrections to lensing for discrete and continuous mass distribu¬ 
tions in GR 

The method of extracting the 1/6”' corrections to the angular dehection can be extended 
to the case of a countinuos mass distribution in the lens plane, as shown in Fig. 9. For this 
purpose we can consider the case of a distributed lens with a surface density 



(B.l) 


In the Newtonian approximation we have the usual relation 



(B.2) 


or, after the rescaling ^ = ^qx, 



(B.3) 


The dehection angle is 



(B.4) 
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and is the gradient of the lensing potential 


a{x) = Vx^{x) 


= — [ (fx' log \x — x'\ 

TT / ^r.r 


where we have introduced the critical surface density 

1 Dos 


^r.r — 


dvr G M Dls Dql 


(B.5) 


(B.6) 


It is possible to introduce the corrections to this behaviour by extending to the case of a 
continuous distribution the result given in Eq. (5.3). Inserting the generalized 1/6 expansion 
given in (5.3), we derive a generalized version of (B.2) that remains valid also in the case 
of strong lensing and which is given by 
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This allows to define a generalized deflection potential 

a{x) = Vx’J'(f), 

with 
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C Full cross section 

We give the explicit expression of the neutrino differential cross section in terms of E and 
9 used for the calculation of the form factors 
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D Coefficients of the semiclassical expansion 

The coefficients present in the expansion in Eq. (4.12) are given by 
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We find 
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E Scalar integrals 

In this appendix we collect the definitions of the scalar integrals appearing in the computa¬ 
tion of the matrix element. One-, two- and three- point functions are denoted respectively 
as ^0) Eq and Cq with 
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iTT^ J [E - mQ){[l +piy - +P2) - "^ 2 ) 

Because the kinematic invariants on the external states of our computation are fixed, = 
{pi — P 2 )^) Pi = P 2 — 0) have defined the shorter notation for the three-point scalar 
integrals 


Co{mQ, m\,m\) = Co(m^, ,m?, ttiq, m\,m \), (E-2) 

with the first three variables omitted. 


F Deflection integral 

The deflection integral for the GR solution can be recast in the form 


0.3/2 _ dt _ 

^ ° Jo 7x3^3 _ ^3^2 + _ 1 


- vr. 


(F.l) 
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with xq being the point of closest approach between the deflector and the beam. It is an 
elliptic integral of first kind. In the indefinite form, the general expression of these types of 
integrals is given by 


In = 


dz 




with the polynomial at the denominator 

R{z) = cqz” + ClH-h c„ 


(F.2) 


(F.3) 


being of degree n. For n = 3, which is the GR case, the Weierstrass form of (F.2) is obtained 
by introducing the roots (ci, 62 , 63 ) and given by 


In = 


dt 


By the transformation 


\/4(f - ei){t - e 2 ){t - 63 ) 


, Cl - 63 

t — 63 -s— 

z^ 


it can be brought to the Legendre form 

Io= ^ 


CdC 


v/(l-C2)(l-fe2C2)’ 

where k is its modulus. In the case of a finite integration, the form that is needed is 

dz 


F{x;k)= r 
Jo 


also re-expressed as 


F{(j^;k)= [ 
Jo 


Y^(l — z‘^){l — k'^z'^) ’ 
d^ 


'0 -s/l — k"^ sin^ d 

by a simple change of the integration variable. For (F.l) the corresponding roots are 


(F.4) 

(F.5) 

(F. 6 ) 

(F.7) 

(F. 8 ) 


ei = — ( xo - 1 - •\/ x5 -h 2x0 - 3 


62 = — ( xo - 1 + yxg -h 2X0 - 3 


63 = — 
Xo 


(F.9) 


and the transformation (F.5) is given by 


^ _ XQ - 3 - v^Xq -F 2x0 - 3 

Xo 2 xo2^ 


which takes to the Legendre form 

Io = -M 


8 x 0 


dz 


Xo 


- 3 - i/xg 2x0 - 3 7 v^(l - z‘^){l - k‘^{xn)z‘^y 


(F.IO) 


(F.ll) 
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with 


(F.12) 


fe^(xo) = 


Xo 


— 3 + 'x/ Xq + 2xo — 3 


Xo - 3 - \/xl + 2xo - 3 
Keeping into acconnt the finite integration region, (F.l) becomes 

dz 


a(xo) = -4S(xo) 


'a(xo) ^(1 - z^){l - hP‘{xo)z‘^y 


with 


A(xo) = 


'Xo 


- 3 - y^Xq + 2xo - 3 


2 x 0 


(F.13) 


(F.14) 


The bonndaries of integration in (F.13) can be expressed in the Jacobi form by the snbsti- 
tntion ^ C = obtaining 

/■°° dz 1 (iQ 

Jam 70^(1 - fe2(xo)z2)" Jo 

with r(xo) = (A(x0))“^ e A^(xo) = (A:^(xo))~^- After reabsorbing a factor in front of the 
integral in the definition of S(xo), we obtain 

a(xo) = -4S(xo)T((/)(xo); A(xo)) - vr, 


with 


S(xo) = 


1 


X 


(^Xo - 3 - Y^xg + 2X0 - 3^ 
6 — 4x0 


4>{xo) = arcsin(r(xo)) 


A^(xo) = 


I Xo - 3 - Y^x§ + 2xo - 3 
Xo - 3 - \/xl + 2xo - 3 


Xo 


— 3 + x/X q T 2xo — 3 


(F.16) 

(F.17a) 

(F.17b) 

(F.17c) 

(F.17d) 


G Feynman Rules 


We collect here all the Feynman rnles involving an external gravitational field that have 
been nsed in this work. All the momenta are incoming 

• graviton - gange boson - gange boson vertex 


V'’ 


= {h • k2 + M^) k2) + fes) 

(G.l) 



where V stands for the vector gange bosons g, 7 , Z and W. 
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• graviton - fermion - fermion vertex 





= It'" (^1 + k 2 Y + i' {ki + k 2 Y - 2 t'"'" (^1 + ^2 - 2m/) I 

(G.2) 


• graviton - scalar - scalar vertex 


, - i|2x|(fci + k2Y{ki + k2Y - V^Yki + ^2)'} 

(G.3) 

where S stands for the Higgs H and the Goldstones 4> ■ The first line is the 

contribution coming from the minimal energy-momentum tensor while the second is 
due to the improvement term. 


• graviton - scalar - fermion - fermion vertex 



(G.4) 


where the coefficients are defined as 


cY.. = cl.:. = -i 


m 


hipxfj 


hipi/j 




2sw Tn\Y 

m^ 


m^ 


Vj 




nL_ _ — A 


m 


.R 


'2swmw 

6 Tfl-ijj 

■ Vj,ni 


2 / 3 , 


^0+bb ^Y2sw^w''^^' 

^ V2swmw 


(G.5) 


• graviton - gauge boson - fermion - fermion vertex 



(G.6) 
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with 


^g'lpip ~ ~ 




C^j , =i - (v + a), C^j I =i - (v — a), 

2 SV1/CVF 2swcw 

e 


a 


w+i^i^ - ^ ^ 0 

and V = I3- 2 s‘^Q, a = I3. 

The tensor structures C, D and E which appear in the Feynman rules defined above are 
given by 


Cfj,iipcr — Vgp Vl'cr T Vpcr 'Hl'p Vpi' Vpcr ) 

D pypa (^1 ) ^2 ) — ^pu klcF ^2 p 1]^ ki /C 2 T Vpp k\a k2u Vpcr k\ ph2u T (M ^ 


(G.8) 


Epvpa (^1 ) ^ 2 ) — Vpv {kl pk\cr T ^2 p ^2 (T T k\ p k2 a) 


VucT kipkip + rj^p k2pk2a + (^ o jv) 
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